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摘要 : 年 轻 脉冲 星 PSR J0614+2229 是 少数 模式 变换 时 发 生 相 位 偏 移 的 脉冲 星之 一 。 本 文 利 
用 Parkes 射电 天 文 台 的 10 cm、20 cm 和 50 cm 三 个 波段 上 的 偏振 观测 数据 ， 对 每 个 波段 爆 
发 模式 和 正常 模式 的 线 偏 振 位 置 角 用 旋转 矢量 模型 进行 拟 合 ,确定 了 其 最 大 斜率 点 的 相位 差 。 
进而 利用 广泛 使 用 的 确定 辐射 高 度 的 相对 论 性 方法 ,首次 得 到 了 模式 变换 时 辐射 高 度 的 变化 
发 现 正常 模式 的 辐射 位 置 高 于 爆发 模式 ,其 高 度 差 随 着 频率 增加 而 扩大 。 根 据 两 个 模式 累积 
脉冲 窗口 中 心 的 相位 差 ， 可 以 推断 两 个 模式 的 辐射 区 并 不 都 关于 磁 轴 -自转 轴 所 在 子午 面 对 


称 。 
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PSR J0614+2229 ( BI B0611+22 ) 的 自转 周期 为 0.335s， 特征 年 龄 只 有 9 万 年 ， 非 常年 
轻 ， 是 一 颗 特 殊 的 模式 变换 脉冲 星 。 除 了 像 其 它 一 些 模式 变换 脉冲 星 那样 有 爆发 Burst , 
以 下 简称 B 模式 ) 和 正常 ( Normal , 以 下 简称 N 模式 ) 两 种 不 同 的 辐射 状态 之 外 ， 其 脉冲 
相位 也 存在 着 同步 的 偏 移 ， 并且 其 模式 变换 呈现 准 周期 性 。 

自 Davies SA"! 于 1972 年 发 现 这 颗 星 以 来 ， 其 模式 变换 行为 已 被 多 个 团队 研究 过 。 
1980 年 ,Ferguson 和 Boriakoff "! 最 早 通过 单 脉 冲 观测 注意 到 其 辐射 存在 时 标 为 几 分 钟 的 变 


化 , Œ 430 MHz 上 ,他 们 发 现 当 这 颗 星 的 脉冲 较 强 时 ， 其 相位 常常 偏 后 一 些 ， 脉冲 轮廓 的 


相位 偏 移 量 可 达 2°72 A. 1992 年 NowakowskiQ 通过 在 相同 波段 研究 ， 确 定 它 为 模式 变换 
脉冲 星 ， 模 式 特征 时 标 为 几 百 秒 。 两 种 模式 下 的 脉冲 轮廓 都 为 高 斯 形 ， 但 宽度 、 强 度 和 相位 
均 有 差别 。2000 年 Ken? 发 现 这 颗 星 的 偏振 位 置 角 (PPA ) 的 最 大 斜率 处 在 脉冲 轮廓 的 尾 
部 边界 ,由 此 认为 其 辐射 主要 来 自 脉冲 星 极 冠 区 ( Polar Cap ) 前 部 区 域 ,而 不 是 充满 极 冠 区 
的 ; 并 提出 B 模式 脉冲 轮廓 更 宽 是 由 于 存在 暂 现 锥 成 分 的 猜测 。2014 年 Seymour SA"! 通 
过 327 MHz 和 1400 MHz 两 个 波段 上 在 同一 天 进行 的 单 脉 冲 观测 发 现 ， 在 327 MHz 上 相位 
较 后 的 模式 明显 强 于 相位 较 前 的 模式 ， 而 在 1400 MHz 上 两 种 模式 的 强度 差别 不 大 。 此 外 ， 
他 们 还 通过 谱 分 析 试 图 寻找 子 脉 冲 漂移 特征 , 但 两 种 模式 都 没有 发 现 子 脉 冲 漂移 的 证 据 。 
2016 Ẹ , Rajwade $A"! 利用 在 150 MHz 和 327 MHz 上 的 同时 观测 ， 以 及 在 327 MHz 和 
820 MHz 上 的 同时 观测 ,发 现在 150 MHz 和 327 MHz 上 的 B 模式 的 相位 都 靠 后 些 ， 而 在 
820 MHz 上 则 相反 ，N 模式 反而 更 靠 后 ， 认 为 Kern 提出 的 暂 现 锥 成 分 假说 尽管 可 以 解释 B 
模式 脉冲 宽度 的 增加 ， 但 无 法 解释 在 不 同 频率 上 B 模式 相对 N 模式 的 相位 变化 。 

以 上 研究 发 现 了 这 颗 脉冲 星 复杂 的 模式 变换 行为 ,但 对 其 模式 变换 时 辐射 区 域 如 何 变化 
还 没 给 出 很 好 的 解释 ,这 很 大 程度 上 是 由 于 过 去 缺乏 多 波段 同时 偏振 观测 。 已 有 的 偏振 观测 
表明 该 星 的 PPA 是 一 条 光滑 的 S 形 曲线 59 ,可 以 用 经 典 旋转 矢量 模型 (RVM rotating vector 
model) “1 或 相对 论 性 RVM 模型 ( RRVM ) "很 好 地 描述 ， 并 以 此 估计 出 了 辐射 高 度 。 
但 是 文献 中 尚 无 对 不 同 模 式 的 偏振 位 置 角 和 辐射 高 度 进行 研究 的 工作 ,这 正 是 本 文 研究 的 重 


点 
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本 文通 过 分 析 Parkes AMG 64 米 射 电 望 远 镜 获 得 的 10 cm 和 50 cm 波段 同时 观测 数据 
以 及 20 cm 波段 的 观测 数据 5 ,从 多 个 波段 考察 其 模式 变换 时 的 偏振 特征 的 变化 ,并 基于 
RVM 模型 对 PSR J0614+2229 辐射 区 域 的 分 布 和 随 频 率 的 变化 给 出 新 的 结果 。 本 文 第 1 部 分 
介绍 了 数据 来 源 和 处 理 流程 ,第 2 部 分 研究 了 两 种 模式 的 辐射 高 度 的 差异 及 其 随 频 率 变化 的 


关系 ， 第 3 部 分 为 总 结 和 讨论 。 


1. 观测 与 数据 处 理 


本 文 使 用 了 Parkes 天 文 台 释 放出 来 的 64 米 射电 望远镜 的 历史 观测 数据 "93。10 om, 
20 cm 和 50 cm 三 个 波段 的 中 心 频率 分 别 为 686 MHz, 1369 MHz 和 3100 MHz , 所 用 后 端 均 
为 WBCORR 采用 1024 采样 数 和 1 min 子 积分 时 长 。 在 686 MHz 和 3100 MHz 通过 5010CM 
接收 机 在 MID 53658.80926 同时 观测 ,观测 时 长 为 51 min ,有 效 带 宽 分 别 为 64 MHz 和 1024 
MHz ,单个 通道 带宽 分 别 为 0.125 MHz 和 1 MHz。 在 1369 MHz 通过 H-OH 接收 机 在 MJD 
53667.80810 进行 观测 观测 时 长 为 18 min 有 效 带 宽 为 236 MHz 单 通 道 带 宽 为 0.125 MHz. 

脉冲 星 的 原始 观测 数据 为 PSRFITSC1 文件 ,本 文 使 用 了 脉冲 星 数据 处 理 专 用 软件 包 
PSRCHIVE "O 来 处 理 。 首 先 要 把 频 域 两 边 信 噪 比 很 低 的 通道 用 paz 模块 切 掉 ， 并 对 数据 
进行 自动 去 噪 ,用 pazi 模块 手动 消 掉 频 域 和 时 域 的 射电 干扰 ( RFI ) ; 第 二 步 是 利用 射电 参 
考 源 Hydra 和 噪声 二 极 管 的 信号 对 数据 用 pac 模块 进行 流量 定 标 、 偏 振 校 准 和 pem 定 标 。 


最 后 再 对 处 理 后 的 偏振 数据 进行 PPA 的 拟 合 等 分 析 。 
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2. 不 同 模式 辐射 高 度 的 计算 
2.1 计 算 辐 射 高 度 的 方法 


计算 辐射 高 度 的 方法 有 多 种 LH ， 其 中 Blaskiewicz SA” 1991 年 ( 以 下 简称 BCW91 ) 
提出 的 相对 论 性 方法 是 两 种 常用 的 方法 之 一 ， 它 用 RVM 模型 来 拟 合 PPA 数据 ,确定 最 佳 拟 
合 曲线 的 最 大 斜率 点 ( SGP steepest gradient point ) 的 相位 中 , 并 计算 它 与 脉冲 窗口 中 心 相 
位 中 .之 间 的 差 值 ， 用 如 下 近似 关系 来 计算 辐射 高 度 r: 


4r 
Adgcw = 中 o — qc = Re ， (1) 


其 中 Ri 为 脉冲 星 的 光速 圆柱 半径 。 


上 述 相位 差 来源 于 磁 层 共 转 所 引起 的 相对 论 性 效应 一 一 光 行 差 效 应 (aberration) 和 辐 


射 高 度 与 方向 不 同 所 引起 的 光 程 差 效应 (retardation)。 相 比 之 下 ， 另 一 种 常用 方法 一 一 通过 


脉冲 宽度 和 辐射 几何 关系 计算 高 度 的 方法 没有 考虑 这 些 因 素 ,因此 相对 论 性 方法 比 纯 几何 的 
方法 更 加 接近 物理 实际 。 

该 方法 适用 于 r <<Rie 的 情形 ,优点 是 计算 公式 简单 ， 且 不 像 辐射 几何 方法 那样 需要 得 
到 较为 精确 的 磁 倾 角 a ( 磁 轴 和 自转 轴 夹 角 ) 和 碰撞 角 B ( 视线 和 磁 轴 夹 角 )。 而 这 两 个 参数 
对 辐射 几何 方法 是 必须 的 ,虽然 它们 可 以 通过 对 PPA 数据 进行 RVM 模型 进行 拟 合 得 到 ,但 
对 于 像 PSR J0614+2229 这 样 轮廓 较 窄 的 脉冲 星 ,所 得 到 的 a 的 误差 较 大 ( 如 BCW91 WAN 
结果 为 apcwo1 = 30°+ 100°), 因此 用 辐射 几何 方法 得 到 的 高 度 误差 较 大 。 而 相对 论 性 方法 
对 磁 倾 角 的 误差 不 敏感 ， 所 以 本 文选 用 这 种 方法 。 

基于 该 方法 ( 即 (1) 式 ) ,BCW91 计 算 了 PSR J0614+2229 在 1.4 GHz 上 的 辐射 高 度 ， 


结果 为 rgcwo1 = 490 140 km@ ( 对 应 于 相位 差 Agcwoi =742° ) 不 过 这 一 结果 是 针 


© 原文 中 由 于 将 周期 误 作 0.25 s 而 得 到 的 值 为 350 kmo 


对 总 积分 轮廓 的 ， 不 能 回答 不 同 模式 时 辐射 高 度 的 变化 。 本 文 已 经 分 别 得 到 了 三 个 频率 上 B 
模式 和 N 模式 的 偏振 轮廓 ， 原 则 上 可 以 利用 该 方法 计算 各 自 的 辐射 高 度 ， 从 而 检验 是 否 发 
生 了 高 度 变化 ,以 及 随 频 率 的 变化 情况 。 但 是 该 方法 有 一 个 前 提 假 设 ， 即 要 求 辐射 区 域 关 于 
磁 轴 和 自转 轴 所 在 平面 ( 以 下 简称 磁 轴 子午 面 ) 对 称 。 经 过 对 比 两 个 模式 脉冲 窗口 中 心 相位 
以 及 PPA 曲线 SGP 的 相对 移动 ， 可 以 推断 至 少 有 一 个 模式 不 满足 这 一 前 提 假 设 ( 详 见 下 一 
TT), 因此 不 能 简单 的 用 公式 ( 1 ) 来 求解 。 

我 们 注意 到 2008 年 Dyk!” 对 BCW91 的 方法 做 了 深入 研究 ， 并 对 公式 ( 1 ) 中 的 相位 
差 给 出 了 更 为 直观 的 理解 ,6 在 r < 0.1RI。 且 极 冠 区 满足 Gouldreich-Julian 电流 密度 的 情况 
下 ， 作 者 将 BCW91 的 RRVM 模型 改写 为 对 惯性 观测 者 参照 系 (IOF ) 更 加 友好 的 形式 : 


cost sina cos(qops — Qe + A Qr) — sin gcosa 

其 中 视线 与 自转 轴 夹 角 z = a + B ,中 os 为 IOF 参照 系 中 观测 到 的 脉冲 相位 ,qe 为 磁 轴 子午 
面 中 辐射 高 度 为 0 处 的 脉冲 相位 ( 作为 参考 相位 )，A 中 .是 高 度 为 + 处 的 辐射 相对 0 高 度 处 
辐射 的 相位 益 等 价 于 轮廓 窗口 中 心 相位 相对 参考 相位 的 偏 移 量 ) 小 ,为 该 处 辐射 的 PPA ， 
由 为 自转 轴 在 垂直 视线 的 平面 中 的 投影 角 。 由 于 光 行 差 和 光 程 差 效 应 的 联合 贡献 ， 导 致 脉 
冲 窗口 中 心中 .相对 参考 相位 中 f 移 动 了 如 下 相位 差 : 

he T (3) 
其 中 负 号 表示 窗口 中 心 前 移 了 。 要 保证 公式 ( 1 ) 成 立 , 则 要 求 SGP 的 相位 中 相对 于 中 后 移 
了 如 下 的 值 : 

hbo = bo— br~ =» (4) 

i 
也 就 是 说 ,公式 ( 1 ) 中 的 相位 差 可 以 理解 为 脉冲 窗口 中 心 前 移 和 SGP 后 移 共同 贡献 的 结果 。 
上 面 两 个 相位 差 中 ，A 中 .依赖 于 辐射 区 关于 磁 轴 子午 面 对 称 的 假设 ， 而 A 中 0 不 依赖 。 


设想 一 个 简单 的 对 比 ,一 种 情况 是 辐射 区 在 磁 轴 子午 面 两 侧 对 称 分 布 , 另 一 种 情况 是 辐射 区 


只 在 磁 轴 子午 面前 面 一 侧 ( 是 前 者 的 一 半 ), 二 者 辐射 高 度 都 相同 ， 那么 后 者 的 PPA 曲线 就 
是 前 者 PPA 曲线 的 一 半 。 既 然 PPA 形状 重合 ,那么 二 者 的 SGP 就 都 一 样 。 不 过 后 者 轮廓 的 
窗口 中 心 显 然 比 前 者 更 加 靠 前 了 。 由 此 可 见 ，SGP 的 相位 与 参考 相位 的 变化 量 能 准确 的 反 
映 辐射 高 度 ， 对 辐射 区 不 对 称 的 情况 也 适用 。 

本 文 的 主要 任务 是 研究 模式 变换 是 辐射 区 高 度 是 否 发 生变 化 。 由 于 参考 相位 是 共同 的 ， 
所 以 就 可 以 利用 两 种 模式 下 SGP 的 相位 差 来 计算 辐射 高 度 之 差 ， 再 结合 脉冲 窗口 中 心 相 位 
的 差别 来 评估 辐射 区 域 偏离 对 称 的 情况 。 采 用 的 工作 步骤 如 下 。 

(1) 拟 合 各 模式 的 PPA 数据 ,得 到 最 佳 拟 合 曲线 的 SGP 相位 ; 通过 SGP 相位 差 来 

计算 各 频率 上 两 个 模式 的 辐射 高 度 差 。 
(2) ”测量 各 模式 脉冲 窗口 的 中 心 相 位 ,通过 模式 间 窗 口中 心 相位 差 以 及 SGP 相位 与 


窗口 中 心 的 相位 差 来 推断 辐射 区 的 变化 。 


2.2 PPA 的 RVM 拟 合 


原则 上 拟 合 PPA 数据 应 该 使 用 RRVM 模型 ,不 过 在 r < 0.1Ri. 的 情况 下 RRVM 模型 和 
经 典 RVM 模型 所 给 出 的 PPA 曲线 形状 一 致 , 这 可 以 通过 对 上 比 公式 ( 2 ) 和 如 下 经 典 RVM 
模型 "函数 关系 看 出 : 


sin a Sin( 中 obs x do) 
sin cos a — cosisinacos(qops — Po)” 


tan( — Yo) = (5) 


其 中 中 o 和 由 ,是 磁 轴 子午 面 中 磁力 线 上 出 来 的 辐射 相位 与 PPA , 中 与 公式 ( 2 ) 中 的 中 是 等 价 
的 ， 同 时 它 也 是 PPA 曲线 上 最 大 斜率 点 ( SGP ) 的 相位 ,有 : 
dp) _ sina 


(=) = ( 动 ~ sing ’ (6) 


max 


因此 , 为 简便 起 见 ， 本 文 仍 采用 经 典 RVM 模型 拟 合 PPA 数据 。 


另 一 个 需要 处 理 的 技术 问题 是 上 式 中 的 四 个 参量 po。、y。、a 和 B 的 拟 合 结果 依赖 于 数据 
质量。 数据 质量 越 高 ， 拟 合 结果 不 确定 度 越 小 。 此 外 ， 对 于 同一 颗 星 的 多 频 辐 射 和 各 辐射 模 
式 ,尽管 go 和 业 o 可 能 有 变化 , 但 a 和 B 的 值 应 该 保持 不 变 。 基 于 这 两 点 考虑 ， 应 该 选择 数据 
质量 最 高 的 偏振 轮廓 来 拟 合 得 到 a 和 B， 并 将 其 固定 ， 用 于 各 模式 的 PPA 数据 的 拟 合 ， 来 得 
到 中 ,和 小 ,。 拟 合 采用 了 PSRCHIVE 中 的 psrmodel 模块 。 它 是 是 在 IAU 协定 09 下 编写 的 基 
于 经 典 RVM 模型 拟 合 四 个 参数 的 模块 ， 程序 允许 同时 拟 合 四 个 参数 ， 也 允许 限定 部 分 参数 
来 拟 合 其 它 参 数 。 

三 个 波段 的 数据 中 10 cm 和 50 cm 的 观测 时 间 最 长 , 均 为 一 小 时 ,但 脉冲 星 在 50 cm 的 
辐射 要 强 得 多 , 因而 50 cm 数据 的 信 噪 比 ( SNR ) 在 三 个 波段 中 最 高 , 其 PPA 的 误差 最 小 、 
相位 范围 最 大 ,这 对 于 得 到 更 可 靠 的 a 和 B 值 至 关 重 要 。 因 此 ,我 们 首先 选取 了 686 MHz 
总 积分 轮廓 的 PPA 数据 来 拟 合 。 得 到 的 结果 如 图 1 所 示 ， 具 体 参 数 为 : a = 83.8? 78.9° 
Al B = —6.0° + 2.2° 和 BCW91 在 1.4GHz 的 结果 在 IAU 协 定 下 09 agew = 30° + 100° 
和 Bgcwo1 = —6°+15°) 在 误差 范围 内 相 吻 合 。 从 图 1 可 以 看 出 最 佳 拟 合 RVM 曲线 和 
PPA 数据 符合 得 很 好 。 本 文 得 到 的 PPA 的 SGP 大 约 在 轮廓 的 最 右 侧 边界 处 ， 这 和 Kern 
的 结果 一 致 。 接 下 来 就 可 以 将 上 面 得 到 的 c 和 B 值 用 来 拟 合 三 个 波段 不 同 模式 下 的 po 和 


Woo 


1 : 686 MHz 上 PPAR RVM 拟 合 。 顶 部 图 框 中 实 线 为 RVM 的 最 佳 拟 合 曲线 | ERIX 
处 为 SGP (qo, 业 0)。 底 部 图 框 分 别 用 实 线 、 短 划 线 以 及 点 划 线 给 出 了 总 轮廓 I、 线 偏振 L 和 


AIR V, 同时 用 虚线 标 出 了 总 积分 轮廓 的 窗口 中 心 所 在 相位 和 10% 峰 值 流量 水 平 。 


Fig. 1: RVM fit of the PPA data at 686 MHz. The top panel shows the best-fit RVM curve, with 
the SGP (hp, Wo) being denoted by the cross section of the vertical and horizontal dotted lines. 


The bottom panel shows the total intensity (I, solid), linear polarization (L, dashed), circular 
polarization (V, dash-dotted), the central phase of the profiles window determined at the 10% 


level of the pulse peak (vertical dotted line) and the 10% level of the peak intensity (horizontal 


dotted line) of the averaged pulse profile. 


2.3 两 种 模式 的 辐射 高 度 


我 们 通过 肉眼 辨别 剔除 各 波段 数据 中 可 能 包含 两 种 模式 的 子 积 分 ( 通常 是 在 从 一 个 模式 
的 序列 变 到 另 一 个 模式 的 序列 时 过 渡 的 子 积 分 ), 并 将 剩余 的 两 个 模式 的 数据 分 开 并 又 加 得 
到 各 自 的 总 轮廓 ， 然 后 对 两 个 模式 的 PPA 进行 拟 合 ， 结 果 如 图 2 所 示 ， 其 中 SGP 相位 分 别 
AEE ( B 模式 ) 和 红色 ( N 模式 ) 的 虚线 代表 。 从 图 中 可 以 看 出 ， 从 B 模式 变换 为 N 模 
式 时 脉冲 轮廓 和 PPA 的 SGP 相位 都 发 生 了 明显 的 向 右 偏 移 ,并 且 SGP 相位 的 偏 移 量 呈现 出 
随 频 率 逐 步 增 大 的 特征 。 当 把 N 模式 PPA 的 SGP 和 B 模式 的 对 齐 后 ( 黄色 和 绿色 ), 两 者 


的 PPA 数据 和 最 佳 拟 合 曲 线 都 吻合 得 非常 好 ， 这 说 明 磁 场 的 几何 位 形 没有 改变 。 


686 MHz 1369 MHz 3100 MHz 


PPA (deg) 


Flux Density (mJy) 


= oe d = TET zis a 3 Sgr e a es 
2: 两 种 模式 线 偏振 位 置 角 的 RVM 拟 合 。 顶 部 图 框 : 蓝 色 和 红色 分 别 代 表 SGP 对 齐 之 前 
的 B 模式 和 NN 模式 数据 ,绿色 和 橙色 的 点 分 别 代表 SGP 对 齐 后 的 B 和 N MRE. RG 
图 框 : 以 峰值 强度 10% 水 平 确定 的 两 种 模式 的 轮廓 窗口 中 心 。 各 类 点 和 线 型 的 意义 与 图 1 


相同 。 


Fig. 2: The RVM fit for the PPAs of the two modes. Top panel: the blue and red curves show the 
best-fit RVM curves of the B and N modes without aligning the SGPs, while the green and orange 


points and curves represent the those with the SGP alignment. Bottom panel: the centers of 
averaged pulse profiles’ window determined at the 10% level of the pulse peak. The meanings of 


the data points and the line types are the same as those in Fig. 1. 


表 1 给 出 了 每 个 波段 上 两 个 模式 bo 和 中 o 的 最 佳 拟 合 值 . 它们 之 间 的 相位 差 sbo = 中 op 一 
don = A Pos — ^ 中 oN 和 高 度 差 6r = re 一 ry ( 二 者 的 关系 遵从 公式 ( 4 ) ) ,以 及 轮廓 窗口 中 
心 的 相位 中 .和 两 个 模式 之 间 的 窗口 中 心 相位 差 6q。 = bop 一 中 cn = A 中 rp 一 A 中 ew， 其 中 下 
标 B 和 NN 分 别 代 表 对 应 的 模式 。 轮 廓 窗口 中 心 相 位 测量 的 方法 是 通过 取 每 个 轮廓 两 侧 强 度 


等 于 轮廓 峰值 强度 10% 所 在 相位 ( 如 图 2 所 示 ) 为 脉冲 窗口 的 边界 8 ， 记 为 Ll。 和 Rio, 


© 10% 峰 值 宽度 能 够 涵盖 绝 大 部 分 的 脉冲 相位 , 相 比 50% 峰 值 宽度 更 能 代表 脉冲 窗口 。 若 采用 50% 宽 度 ， 
则 反映 的 是 辐射 区 中 较 强 的 区 域 ， 但 是 受到 脉冲 轮廓 形状 的 影响 ,50% 强度 的 脉冲 窗口 和 10% 的 窗口 未 必 
同时 关于 磁 轴 子午 面 对 称 。 


脉冲 窗口 的 中 间 相 位 为 中 .= (Lio 十 Ri0)/2 ; 单 侧 相位 的 误差 可 以 用 公式 a( 中 age) = 
o(Iom)/1d1/d 中 |eage 得 到 ， 其 中 o(lor) 代 表 噪声 区 域 强度 的 标准 方差 ， 导 出 量 的 误差 通过 误差 


传递 公式 计算 。 计 算 中 用 到 了 脉冲 星 的 光速 圆柱 半径 Ri. = 1.6 x 104 kmo 


表 1 : 两 种 模式 的 相位 差 与 高 度 差 


Table 1: The phase offset and the difference in altitude between the two modes. 


vMHz Mode bo/。 Wol? 6q0/® or/km A? 5,/° 
B 90.7(2)  39.3(9) 83.3(3) 

686 f -0.5(3)  -70(40) -0.6(4) 
N 91.2(2) 41.19) 83.9(2) 
B 284.1(2)  -30.8(9) 276.8(3) 

1369 -0.8(5)  -110(70) -0.9(6) 
N 284.9(5)  -28(2) 277.1(5) 
B 299.15)  57(2) 293.8(4) 

3100 -2(1) -270(150) -0.8(6) 
N 301(1) 64(6) 294.6(4) 


注 : 括号 中 为 参数 末 位 的 误差 值 。 
Note: The values in the brackets stand for the errors in last digits. 

表 中 负 的 6qo 值 表示 N 模式 的 SGP 相位 比 B 模式 的 更 靠 后 ， 对 应 于 负 的 Sr 值 ,表示 N 
模式 辐射 高 度 比 B 模式 的 高 。 若 两 个 模式 的 辐射 区 都 对 称 ， 由 于 公式 (3 ) 和 (4) 反 号 ， 
则 8 中 。 应 该 为 正 值 ,这 就 要 求 B 模式 的 轮廓 窗口 中 心 要 比 N 模式 的 滞后 ,但 观测 结果 恰好 相 
反 ,这 说 明 ( 3 ) 式 成 立 的 前 提 一 一 辐射 区 关于 磁 轴 子午 面 对 称 不 可 能 对 两 个 模式 同时 成 立 。 

为 了 解释 观测 N 模式 轮廓 窗口 中 心 滞后 的 现象 ， 理 论 上 存在 三 种 可 能 的 组 合 ,( 1 ) 如 
果 B 模式 的 辐射 区 是 对 称 的 ， 那 么 N 模式 的 辐射 区 就 不 对 称 ， 且 位 于 磁 轴 子午 面 后 侧 的 部 
分 要 比 前 侧 更 宽 一 点 ( 2 ) 如 果 N 模式 的 辐射 区 是 对 称 的 ,那么 B 模式 的 辐射 区 就 不 对 称 ， 
且 在 磁 轴 子午 面前 侧 的 部 分 比 位 于 后 侧 的 部 分 更 宽 ;( 3 ) 两 种 模式 的 辐射 区 都 不 对 称 ,但 是 
要 求 B 模式 的 辐射 区 要 上 比 N 模式 更 靠 前 。 

第 一 种 和 第 二 种 情况 下 至 少 有 一 个 模式 的 辐射 区 被 假设 为 对 称 的 ,这 可 以 方便 地 计算 高 


度 。 对 第 一 种 情况 ， 我 们 测量 了 三 个 波段 下 B 模式 的 SGP 和 轮廓 窗口 中 心 的 相位 差 


Adscws = Dog 一 Pcp ( 等 式 右 边 的 两 个 相位 值 见 表 1), 并 用 公式 (1) 计算 了 辐射 高 度 ， 
然后 利用 表 1 中 的 高 度 差 计算 了 N 模式 的 辐射 高 度 。 用 同样 的 办 法 也 计算 了 第 二 种 情况 下 
的 两 种 模式 的 高 度 。 结 果 如 表 2 所 示 。 综 合 表 1 和 表 2 的 结果 ,可 以 看 出 两 个 模式 之 间 的 高 
度 差 随 频率 递增 ,但 是 具体 到 每 个 模式 的 高 度 随 频率 的 变化 关系 , 却 存在 不 确定 性 : 在 假设 
B 模式 的 辐射 区 是 对 称 的 情况 下 ,B 模式 的 辐射 高 度 随 频率 递减 ,而 N 模式 却 是 微弱 递增 ; 
若 假设 N 模式 的 辐射 区 是 对 称 的 ， 则 两 种 模式 的 辐射 高 度 都 随 频率 递减 。 

最 后 , 为 了 方便 对 比 ， 我 们 也 仿照 BCW91 计算 了 各 频率 下 累积 脉冲 总 轮廓 ( 两 种 模式 
混合 FY ied ie AH SGP 和 轮廓 窗口 中 心 相 位 差 累 积 脉冲 轮廓 Apcwr 和 辐射 高 度 ri 如 
表 2 所 示 b 从 结果 看 ， 总 积分 轮廓 得 到 的 辐射 高 度 随 频率 呈现 递减 的 趋势 ， 高 度 值 更 接近 


被 假设 为 辐射 区 对 称 的 那 种 模式 ( 第 一 种 情况 的 B 模式 或 第 二 种 情况 下 的 N 模式 b 


表 2 : 脉冲 星 PSR J0614+2229 不 同 模式 的 辐射 高 度 . 
Table 2: The emission heights of PSR J0614+2229 for different modes. 


v/MHz 4 中 bcwB/” TB 1/ km TN1/ km Adsgcwn/° TN 2/ km TB 2/ km Adgcwr! ° FT/ km 


686 7.4(4) 520(30)  590(50) ”7.3(3) 510(20)  440(50) 7.5(1.0) 520(70) 
1369 7.3(4) 510(30)  620(80) ”7.2(7) 500(50)  390(90) “7.4(2) 520(20) 
3100 5.3(6) 370(40)  640(160) 6.4(1.1) 450(80) 180(170) 5.6(5) 390(30) 


注 : 括号 中 为 参数 未 位 的 误差 值 。 


Note: The values in the brackets stand for the errors in last digits. 


3. 结论 


本 文 利用 2005 F Parkes 天 文 台 在 10cm、20cm 和 50cm 三 个 波段 的 偏振 观测 数据 ， 
对 年 轻 脉冲 星 PSR J0614+2229 模式 变换 时 的 辐射 区 磁场 结构 和 辐射 高 度 可 能 的 变化 进行 
了 研究 。 通 过 使 用 RVM 模型 对 B 模式 和 N 模式 的 偏振 位 置 角 数 据 进 行 拟 合 ,发现 两 种 模式 


的 PPA 数据 和 最 佳 拟 合 模型 曲线 都 很 好 的 吻合 ， 表 明 在 模式 变换 时 脉冲 星 磁场 结构 保持 不 


变 。 根 据 确 定 辐射 高 度 的 相对 论 性 方法 以 及 三 个 波段 上 N 模式 PPA 曲线 的 SGP 相位 相对 于 


B 模式 的 滞后 量 ， 计 算出 了 两 种 模式 高 度 之 差 在 70+40 km 到 270+150 km Zi , BRM 


率 增加 而 递增 。 在 相对 论 性 RVM 模型 框架 下 ， 测 量 得 到 的 两 种 模式 脉冲 窗口 的 中 心 相位 之 


差 并 不 能 与 它们 的 SGP 相位 差 相 自治 ,因此 可 以 推断 这 两 个 模式 的 辐射 区 并 不 都 关于 磁 轴 


子午 面 对 称 分 布 。 但 受 限 于 数据 的 信 噪 比 , 本文 没有 对 脉冲 窗口 变化 对 应 的 辐射 区 经 度 范围 


的 变化 进行 讨论 , 有 待 进一步 的 偏振 观测 获得 更 高 质量 的 多 频 ( 或 超 宽 带 ) 观测 数据 。 
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On the Multi-frequency Emission Regions of the Mode 


Switching Pulsar PSR J0614+2229 
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2. Yuncheng University, Yuncheng 044000, China; 3. Central China Normal University, Wuhan 430079, China.) 


Abstract: The young pulsar PSR J0614+2229 is one of a few known pulsars with phase offset 
during mode switching. With the polarimetry data at 10 cm, 20 cm and 50cm archived at Parkes 
Observatory, we employed the Rotating Vector Model to fit the linear polarization position angle 
data for the burst and normal modes at each frequency and obtained the phases offsets between the 
steepest slope-rate points of best-fit curve. Using the relativistic method that is widely used to 
determine the emission altitude, we firstly calculated the difference between the emission altitudes 
of the two modes. It was found that the emission altitude of the normal mode is higher than that of 
the burst mode, and the difference increases with frequency. On the basis of the phase offset 
between the centers of the averaged pulse windows of the two modes, it was inferred that the 
emission regions of the two modes should not be symmetric, in the same time, about the meridian 


plane that contains the rotation and magnetic axes of the pulsar. 
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